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Abstract 

Using a nonperturbative weak noise approach we investigate the interference of noise and chaos 
in simple ID maps. We replace the noise-driven ID map by an area-preserving 2D map modelling 
the Poincare sections of a conserved dynamical system with unbounded energy manifolds. We 
analyze the properties of the 2D map and draw conclusions concerning the interference of noise on 
the nonlinear time evolution. We apply this technique to the standard period-doubling sequence 
in the logistic map. From the 2D area-preserving analogue we, in addition to the usual period- 
doubling sequence, obtain a series of period doubled cycles which are elliptic in nature. These 
cycles are spinning off the real axis at parameters values corresponding to the standard period 
doubling events. 
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I. INTRODUCTION 



The noise-induced escape from an equilibrium state is a fundamental problem in many 

areas in science. The classical work goes back to Kramers ll|l who computed the transition 

□ 

rate from a single potential well; see also the review in Ref. More recently, the under- 
standing and description of nonlinear dynamical systems exhibiting e.g. period-doubling, 
intermittency, and chaos, have led to a renewed interest in the influence of noise on the 

In this context the generation of chaotic motion in low dimensional time-discrete systems 
P.[lo[ 11 , 1^ offers a particularly simple case lending itself to the analysis of the influence of 
noise. Thus a variety of phenomena have been investigated such as the noise-induced shift, 
broadenin g, a nd suppression of bifurcations 1^ 
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and the invariant density |17l . 
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br chaos, the scaling 
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see also recent work on the application of periodic orbit theory to noisy maps 

For the purpose of studying the interference of stochastic noise with the nonlinear be- 
havior of dynamical systems, the simplest case is the influence of noise on 1-D maps of the 
generic type 



n+l 



Here n is the discrete time index and the map f{x) defines the nonlinear discrete time 
evolution. In the linear case /(x) oc x the map is readily analyzed; it possesses an attractive 
or repulsive fixed point at x = depending on the slope, and does not exhibit chaos. In the 
nonlinear case, where /(x) possesses one or several differentiable maxima, the map generates 
the well-known period-doubling sequence to a chaotic state beyond a critical value of the 



control parameter characterizing the map 
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. In the singular case of for example the 



tent or the shift map, the period-doubling sequence is absent and the time evolution becomes 
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chaotic beyond a critical control parameter value, see e.g. Refs. 

Assuming that the deterministic map in Eq. (jl.ll) describes the dynamical evolution of 
a nonlinear system in an effectively reduced phase space, the presence of stochastic degrees 
of freedom originating from the environment can typically be modelled by an additive noise 



2 



term and we are led to the stochastic map 



f{Xn) + in\ 



;i-2) 



a discrete version of the usual Langevin equation for stochastic processes 

Here ^„ indicates the applied noise which we for simplicity assume has a Gaussian distri- 
bution 



V2^ 



exp 



Sn 

2A 



;i.3) 



and is correlated according to 



:i.4) 



where A is the noise strength and (■ ■ ■) indicates an average over the noise ensemble. 

The map / acts like a nonlinear filter on the known white noise input fluctuations C,n 
and the fundamental issue becomes that of determining the stochastic properties of the 
stochastic variable x„. The natural small parameter is the noise strength A which enters 
in a nonperturbative manner as indicated by the form of -P(0 Eq. ()1.3|) . Physically, the 
case of vanishing noise A = 0, yielding the dissipative map in Eq. (jl.lj) for ,^ = 0, is quite 
distinct from the case of weak noise A ~ 0, where on a sufficiently long time scale (labelled 
by the iteration index n) the noise ^„ drives the variable x„ into a stochastic state. The 
issue is to understand how the stochastic noise interferes with the nonlinear behavior. 

In recent work we have elaborated on a nonperturbative canonical phase space approach 
based on the Freidlin-Wentzel formulation js^ or, alternatively, a saddle point approximation 
to the Martin-Siggia-Rose method in its functional form js^ls^ and have applied this scheme 



to the stochastic 
interface 
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K^ardar-Parisi- Zhang equation jail 
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l38l | describing the evolution of a growing 



It is characteristic of the canonical phase space formulation that the stochastic evolution 
equation is replaced by coupled Hamiltonian equations of motion, the noise being replaced 
by a canonical momentum variable, and that the transition probability is obtained from 
the action associated with a solution (an orbit in phase space) of the Hamilton equations. 
The formulation is a weak noise approximation in the same spirit as the well-known WKB 
approximation in quantum mechanics and the noise strength enters in the same way as the 



Planck constant in the WKB scheme 



45|. 



In the present paper we attempt to adapt the above singular weak noise scheme to a 
stochastic map. We again find a characteristic variable doubling in the sense that the 
stochastic map in Eq. (|1.2|) is replaced by the 2D map 



f{Xn) +Pn, 
Pn~l = f'{Xn)Pn, 



:i.5) 
:i.6) 



where the noise C,n in Eq- ()1-2|) is replaced by the new noise variable pn- Note also that the 
equation for p„ iterates 'backwards' in time. The transition probability from Xn to Xm in 
time A^, P(x„ — > Xm, N) is given by 



with dynamic partition function 

n{N) = J^exp 

Xm 

The action S has the form 



:i.7) 



A 



1 ^ 



[1.9) 



n=l 



The formal scheme is thus straightforward. From the 2D map given by Eqs. and ()l.fi|l 
we extract an orbit from a;„ to Xm traversed in steps. The initial and final condition on 
X together with the time span then defines Pn and P follows from Eq. ()1.7|) by means of 
the action 5* evaluated in Eq. ()1.9p . 

II. WEAK NOISE SCHEME 



A. Generalities 

Generally, the stochastic properties of the noisy map in Eq. (jl.2p can be extracted from 
the generator 



Y\_dxr. 



(2.1) 
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where Xn is driven by the map and (■ ■ ■) indicates an average over the imphcit noise depen- 
dence. Implementing the map by means of a deha function constraint which is subsequently 
exponentiated, and averaging over the noise according to Eq. (|1.3p we obtain 

•' n 

Finally, integrating over the noise variable ^„ and making the replacement iAp„ — + — p„ we 
arrive at the form 

Z({/i„}) OC I J]rfx„rfp„e'^"^"e-^/^, (2.3) 
where the action S is given by 

S = ^Pnixn+l- f{Xn) -]^P^ . (2.4) 

In the weak noise limit A ^ the dominant part of the functional integral in Eq. (|2.3p 
is determined by the orbits minimizing the action 5* and we obtain a principle of least 
action yielding the condition 5S = Q subject to variations of x„ and p„. Setting 6S/6xn = 
Pn-i — Pnf'{xn) = and 6S/6pn = Xn+i — f{xn) — ]5n = wc then obtain the equations of 
motion in the form of the 2D map in Eqs. ()1.5|) and p.6p . For the optimal orbit traversed 
in time N we also obtain by inserting Eq. ()1.5p in Eq. ()2.4j) the action in Eq. p.9|) . From 
Z{{fin}) we extract the transition probability from x„ to Xm in time N given by Eq. (jl.7p . 

In order to evaluate the transition probability from Xn to Xm in time N we then have 
to first solve the 2D map given by Eqs. ()1.5|) and ()1.6p subject to the conditions Xq = Xn 
for g = 1 and Xq = Xm for q = N, pq being a slaved variable; secondly, evaluate the action 
S for this specific orbit, and finally evaluate P according to Eq. ()1.7|) . We note that the 
delta function constraint implementing the map gives rise to the additional noise variable 
Pn- We also remark that the equation for pn iterates backwards in time. Solving for pn we 
can express the 2D map in a form iterating forward in time 

Xn+1 = f{Xn) +Pn, (2.5) 
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Eliminating p„ from Eqs. ()2.5|) . ()2.6p . and inserting in Eq. p.9|) we can also express 5* in 
the form 



1 ^ 



(2.7) 



n=l 



in combination with the two-step recursion formula 



~ fiXn-l) ~ i^n+l ~ /(^n))/ i^n) — 0, 



(2.8) 



thus making contact with the work in Refs. |l3l. Il4l Il7|. 



B. Continuum limit - Symplectic structure 



It is instructive to perform the continuum limit x„ — > x(t), Pn — > p(t). Setting x„+i ~ 
X + dx/dt and Pn-i ~ p — dp/ dt we obtain from Eqs. p.Sj) and p.6|) the coupled differential 
equation, the 2D flow 

dx 



dt 



-X + f{x) +p, 



(2.9) 
(2.10) 



originating from the Hamiltonian 



H = p{-x + f{x)) + -p^; 



(2.11) 



we note that compared with ordinary mechanics the Hamiltonian has a momentum depen- 
dent potential and unbounded energy surfaces in phase space. Correspondingly, the action 
in Eq. ()1.9|) takes the symplectic form 



S 



dt 



dx ■ 



(2.12) 



We conclude that in the continuum limit the weak noise approach to the stochastic flow or 
Langevin equation 

dx 



dt 



-x + fix)+^, 



(2.13) 



with noise correlations 



{mmm = A5it-t'), 
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(2.14) 



yields a canonical structure described by the conserved 2D flow in Eqs. ()2.9p and ()2.10p . 
The phase space plot is spanned by x and p and the orbits lie on the energy surfaces given 
by Eq. (j2.1H) . The zero-energy surfaces are spanned by p = 0, the transient manifold, and 
p = 2x — 2f{x), the stationary manifold. The transition probability from Xi to X2 in time 
T is given by P oc exp[—S/A] and is evaluated by solving the equations of motion (|2.9p 
and ()2.10|) for an orbit from Xi to X2 traversed in time T and computing the action S in 
Eq. ()2.12|) . At long times the orbit migrates to the zero-energy manifolds and pass by the 
saddle point at the intersection of the two manifolds yielding an ergodic stationary state. 
Expressing the Langevin equation in the form 



dx _ IdF 

dt 2 dx ' 



(2.15) 



where the free energy is given by 



X 



f{y)dy, 



(2.16) 



it follows that the system can attain a stationary state with probability distribution 



Pst oc exp 



F 

'a 



(2.17) 



This result also follows from S = J dtpdx/dt on the zero-energy manifold setting 
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2x—2f{x). We shall not pursue this analysis further here but refer to Refs. 
for more details. The main point is that in the continuum limit the weak noise approach 
applied to a noisy ID flow yields a symplectic structure with conserved 2D flow. On the 
other hand, the 2D flow does not exhibit chaotic behavior which is the main issue under 
investigation here. 



C. General discussion of 2D map 

The forward iterating 2D map in Eqs. ()2.5|1 and ()2.(jj) is not area-preserving. This 
follows from the Jacobian J = d{xn+i,Pn+i)/d{xn,Pn) = fi^n)/ f'{xn+i) which only in 
the continuum limit approaches unity, implying area-preservation as discussed in Sec. IIIBI 
Only in the vicinity of the flxed points where the orbits slow down does the continuum limit 
apply and we obtain the conservation of area. However, introducing the new noise variable 
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TTn+i = Pn shifted by one itration and coinciding with p„ in the continuum hmit we obtain 
the 2D map 

Xn+l = f{Xn) + -JJJ—^^n, (2.18) 
TTn+l = T77 ^'^n, (2.19) 

which has Jacobian J = 1 and is thus area-preserving. 



1. Fixed Point Structure 

The fixed point structure of the 2D map follows from the equations 

= /(^*) + 7;^' (2-20) 

71* = (2.21) 

with solutions 



(x*,7r*) = (xt,0), (2.22) 
(x*,7r*) = (x;,7r2*) = (x;,x;-/(x;), (2.23) 

where fixed point of the deterministic ID map and solution of /'(xg) = 1 (for 

simplicity we are assuming only one solution), i.e., 

xl = f{xl), (2.24) 
fix;) = 1. (2.25) 

The 2D conserved map is thus characterized by at least two fixed points whose stability we 
proceed to analyze. 

Expanding about a fixed point {x*, vr*) by setting x„ = x* + 6xn and vr^ = vr* + (57r„ we 
arrive at the tangent map j45| 



/ SXn+1 1 




( 5Xn \ 











with stability matrix 



(2.26) 



f'(x*) - TT*/^ 1 
M = r I ^'(-*) I . (2.27) 



Owing to area-preservation the determinant |M| = 1 and the properties of the fixed points 
are determined by the trace 

1 f"(x*) 

fix*) f'{x*y 

and the eigenvalues of M are then given by 




X± = -±\l{-j -1. (2.29) 

The eigenvalues come in reciprocal pairs, A+ = Al^. For |t| < 2 the eigenvalues form complex 
conjugate pairs on the unit circle and the fixed point is elliptic, for t > 2 the fixed point 
is hyperbolic, and for t < —2 the fixed point is inversion hyperbolic. In the limiting case 
|t| = 2 the eigenvalues are degenerate, corresponding to the parabolic case. 

In the case of the deterministic fixed point (x*,7r*) = (x^,0), x^ = f{xl) we have the 
trace, eigenvalues, and invariant manifolds 

t = f'ixl) + (2.30) 

A± = f\x\t\ (2.31) 

((5x,57r) oc (1,0), (2.32) 

((5x,57r) oc (1,1 - /'(x*)2). (2.33) 

Since |t| > 2 the fixed point is hyperbolic or inversion hyperbolic; the degenerate case |t| = 2 
for |/'(xj)| = 1 corresponds to the parabolic case. In Fig. ^ we have depicted the orbits in 
the vicinity of the fixed points in the two cases |/'| > 1 and |/'| < 1. We note that the orbits 
along the invariant deterministic manifold (1,0) corresponds to the ID map. In Fig. |21we 
have shown the trace t as a function of the slope f'{x\) of the ID map. 

For the noisy fixed point {x*,ti*) = {x2,7!'2), vr2 = X2 — f{x2), /'(^2) = obtain, 
setting /2 = /2 (^2)5 trace, eigenvalues, and manifolds 

t = 2- Tilfl (2.34) 

A± = ^ ± (2-35) 

(5x, 5vr) cx 1^1, ^ + + ' (2-36) 

(fe, 57r) oc ( 1, - xlCW^+^n , (2.37) 




FIG. 1: We depict the orbits in the vicinity of the deterministic hyperbohc fixed point (x*,0). In 

a) we show the orbits in the case < 1, where the fixed point for the ID map is stable. In 

b) we show the orbits for > 1; the fixed point for the ID map is unstable. 

and the character of the fixed point depends on the specific map f{x). In Fig. |21we have 
depicted the orbits in the vicinity of the noisy fixed point in the two cases of a hyperbolic 
and elliptic fixed point. In Fig. E] we have shown the trace t as a function of 7r2/"(x2). 



D. Transition probabilities in the vicinity of the fixed points 

Leaving aside for a moment the issue of the nonlinear breakdown of the orbit structure 
near the fixed point the weak noise scheme can be applied in order to determine the transition 
probabilities. 



1. The deterministic fixed point 

In the case of the deterministic fixed point (x*,7r*) = (x*,0), setting a = f'{xi) and 
inserting the initial and final values 6xi and 6xn and the time span A^, the solution of the 
tangent map in Eq. (j2.26|) is 



(2.38) 

6n^ = ia -1)— (2.39) 
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FIG. 2: We plot the trace t of the stabihty matrix M as a function of the slope f'{xl) for the 
deterministic fixed point (x|,0). The shaded area corresponds to the stability region of the ID 
map. 

describing an orbit from 6xi to 6xn in time with SiTn as a slaved variable. 

For \a\ = \f'{xi)\ < 1 the fixed point is stable on the transient invariant manifold 57r„ = 0. 
In the absence of noise the motion on the manifold corresponds to damping, i.e., the ID 
map is line contracting. In the presence of noise for (57r„ ^ the orbits escape from the 
fixed point and approaches the stationary invariant manifold Srcn = (1 ~ CL'^)Sxn', this is the 
scenario depicted in Fig. la). In the continuum limit this behavior corresponds to the flow 
for the noise-driven overdamped oscillator discussed in Ref. Q]. 

For \a\ > 1 the fixed point is unstable on the invariant manifold ^vr^ = 0. However, 
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inverse 
hyperbolic 

FIG. 4: We plot the trace t of the stabihty matrix Af as a function 'ir2f"{x2)- For t > 2 the fixed 
point is hyperbohc, for t < —2 the fixed point in inverse hyperbohc. For |t| < 2 the fixed point in 
elhptic. 

we note that in this case (57r„ 0. The fixed point is attractive along the invariant noisy 
manifold (57r„ = — (a^ — l)6xn] the orbit structure is depicted in Fig. lb). In the vicinity 
of the hyperbolic fixed point the orbit slows down and an orbit from 6xi to 6xn in time N 
must pass asymptotically through the fixed point for oo. 

In order to evaluate the transition probability from 6xi to 6xn in time we insert 
(57r„ = 6pn-i in the expression for the action in Eq. ()1.9j) . Performing the sum yields the 
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action 



D = -(a — 1)(1 — a ^ ' 



7V-1\ 2 



(2.40) 



and from Eqs. p.7|) and p.Sj) the normalized transition probability 

(a^ - 1)(1 - a-2^) /5x;v - feia^-i^ 



P((5xi ^ (5xiv, AT) = fi(A^)"' X exp 
with dynamic partition function 

n(iV) = (a^-i - a-^+i) 



2A 



27rA 



1/2 



(2.41) 



,(a2-l)(l-a-2^). 

In the long time limit N ^ oo the behavior of S, P, and f2 depends on a 
|a| < 1 we obtain 

2 



S 



and thus the stationary distribution 



Pst oc exp 



1 



(2.42) 



/'(x*). For 



(2.43) 



2 A 



5a; AT 



(2.44) 



2. The noisy fixed point 



The noisy fixed point (x*,n*) = (x2, X2 - f (x^)) , fix^) = 1 is elliptic for |7r2/"(x2) - 2| < 
2 and hyperbolic for \'n'2f"i^2) ~ 2| > 2. In the hyperbolic case the discussion in the 
deterministic case applies, we leave the details to the reader. In the elliptic case the complex 
eigenvalues are determined by Eq. ()2.35j) . The orbits are periodic and given by 



6xn = A cos^un + ipi), 
67in = B cos(u;n + ^/'2), 



(2.45) 
(2.46) 



where the frequency uj is given by 



tancu 



2 - vrl/^' 



(2.47) 



In the hyperbolic case, the long time orbits near the fixed point permits a derivation 
of the transition probability as in the case of the deterministic fixed point on the tt = 
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manifold. In the elliptic case, the finite time orbits also allow an evaluation of P; we shall, 
however, not pursue this analysis here. 

Summarizing, the 2D map describing the noisy ID map within the nonperturbative weak 
noise approach is characterized by at least two fixed points. The noiseless deterministic fixed 
point is hyperbolic. The character of the noisy fixed point (s) depends on the magnitude of 
T^'ifi^V)- ^he vicinity of the hyperbolic fixed point the orbits slow down and we can 
derive a stationary probability distribution in the linear regime. Near the elliptic fixed point 
the finite time orbits only yield the short time character of the probability distributions. 

Owing to the nonlinear character of the 2D map we anticipate that the orbit structure near 
the hyperbolic fixed point breaks down due to heteroclinic intersections and the formations 
of foliations. In the vicinity of the elliptic fixed point we expect to encounter the dissolution 
of the KAM tori into cantori and chaotic seas with Arnold diffusion 3G , 45l . Some of these 



features will be investigated in more detail in the next section on the noisy logistic map. 



III. ID NOISY MAPS 



In this section we embark on a discussion of specific noise-driven ID maps which in 
the absence of noise exhibits nonlinear behavior such as period doubling and transitions to 
chaotic behavior. 



A. The linear map 

It is instructive to briefiy consider the linear map 

Xn+l = rXn, (3-1) 

with control parameter r. The Jacobian J = dxn+i/dxn = r and for r < 1 the map is 
length-shrinking with a stable fixed point at x* = 0; for r = the map is length-preserving 
with a line of fixed points. For r > the map is length-expanding, the fixed point at x* is 
unstable and the orbit recedes to infinity. In the presence of noise we obtain the noisy map 

Xn+l = rXn + ^71) (3-2) 

corresponding to the 2D area-preserving map 

Xn+i = rxn + r"V„, (3.3) 
14 



vr„+i = r (3.4) 
This map possesses a hyperbolic fixed point at (x*, vr*) and the discussion in Sec. II applies. 

B. The logistic map 

We now turn to the well-known ID logistic map 

Xn+l = rXn{i — Xn), (3-5) 

which has fixed points at x* = and x* = 1 — 1/r. For r > 3 the fixed point bifurcates 
to period-2 fixed points and in the regime 3 < r < r^, ~ 3.57 the map passes through a 
period-doubling scenario exhibiting universal scaling behavior 111 Il2|. 
In the presence of noise we have the stochastic map 

Xn+l = rXn(l - Xn) + ^n, (3.6) 

and the associated area preserving 2D map 

Xn+l = ''^Xn{l — Xn) H 7" r T, (3-7) 

r(l - 2Xn) 

r(l - 2Xn) 

Since f{x) = rx{l — x), f'{x) = r(l — 2x), and f"{x) = — 2r and according to the general 
discussion in Sec. II the 2D map has two deterministic fixed points on the x- axis and a 
noisy fixed point in the lower half x — tt plane: 

(a;*,7r*) = (0,0), (FPl) (3.9) 

(a;*,7r*) = (l-l/r,0), (FP2) (3.10) 

(x*,7r*) = ((r- l)/2r, -(1 -r)V4r). (FP3) (3.11) 

The trace and eigenvalues characterizing the stability and properties of the fixed points are 
given by Eqs. ()2.28|) and ()2.29|) : likewise, the invariant manifolds follow from {6x,6tt) oc 
(1, — (/'^ — nf" — f'\±)). The fixed point form a triad. The traces associated with the fixed 
points are given by 

ti=r + -, (3.12) 
r 

t2 = 2-r + ^, (3.13) 
t3 = 2-^-^. (3.14) 
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We can now read off the character of the fixed points and establish the correspondence 
with the bifurcations of the noiseless ID map. First, we infer that for all r > (we only 
consider positive r) the deterministic fixed points FPl and FP2 on the x-axis are hyperbolic, 
whereas the noisy fixed point is elliptic for r < and hyperbolic for r > r^, where = 
1 + ~ 3.8284 corresponds to the period 3 tangent bifurcation for the ID map. 

Let us first concentrate at the first period doubling point at r = 3. The Jacobian matrix 
fl2.26|) for the 2D map of the logistic equation is of the form 

r - 2rx* + 



2-K* 1 
2a;*)2 r—2rx* 

For the fixed point FP2 (where tt* = 0) at r = 3 it then takes the form 

M = \ \ . (3.16) 

This corresponds to a degenerate node where the two eigendirections coincide along the 
X-axis. For r < 3 the eigenvalues are Ai < —1 and A2 > —1, respectively. At r = 3 the 
eigenvalues collide in -1 signalling that two period-doublings are about to take place. In 
Fig. El we show this behavior just above r = 3, at r = 3.004. The black dot represents the 
(unstable) fixed point FP2, which seen in the 2D plane now is a hyperbolic fixed point. The 
crosses determines the "usual" stable period two cycle. Simultaneously, a period two cycle 
of elliptic points has been spinning off. Initially, along the eigendirection in the x-axis but 
quickly moving nonlinearly into the plane. We have numerically followed the upper elliptic 
two-cycle point as far as we could for increasing r. The KAM surfaces around it gradually 
dissolves and to the best of our knowledge without the cycle point undergoes further period- 
doubling's. Fig. ini shows this two-cycle point at r = 3.608, with only a few remaining KAM 
surfaces around it, seen on a tiny scale. Furthermore, for increasing values of r, what 
happens is a continuous spinning off of elliptic cycles in close association with the standard 
period-doubling sequence on the x-axis. We show this in Fig. I?] (for the parameter value 
r = 3.462) for the period-doubling to the four-cycle occurring at r ^ 3.449. Note again the 
normal stable four-cycle, represented by crosses, on the x-axis and the associated four-cycle 
of elliptic points in the plane. To substantiate this picture even further, we finally study 
the behavior just above the period-doubling into the stable eight-cycle. At r = 3.55, this 
is shown in Fig. |H1 with the same structure of the stable eight-cycle on the axis with the 
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0.65 0.66 0.67 0.68 0.69 

Xn 

FIG. 5: The phase plane (xn,vr„) for r = 3.004 just above the period-doubhng into the two-cycle. 
The black dot is the now (hyperbolic) unstable fixed point; the crosses are the two-cycle on the 
X-axis and the KAM curves encircle the elliptic two-cycle. 

elliptic eight-cycle in the plane. We thus conclude, that this process will continue in the 
same fashion, with elliptic 2" cycle moving away from the x-axis into the plane at each 
bifurcation point. 

For r > Tc the ID map approaches a strange attractor with positive Liapunov exponent. 
In the 2D map this is reflected by orbits approaching the x-axis subject to infinitely many 
foliations in order to preserve the area. 

IV. SUMMARY AND CONCLUSION 

We have attacked the old problem of noise in ID maps, in particular the logistic equation, 
in a new way. Starting out with a discrete version of the standard Langevin equation, 
we apply a nonperturbative weak noise approach to map the stochastic equation onto a 
deterministic 2D map. This map is area-preserving and the added dimension plays the role 



17 




FIG. 6: One of the elliptic two-cycle points for r = 3.6084 in the "chaotic sea" of the energy 
surface. The KAM surfaces remain only on very small scales close to the cycle point. 

of the noise field. We have derived general properties of this 2D map and applied it to two 
well know systems, the trivial linear map and the logistic map. In both cases we find, that in 
addition to the standard fixed points in the usual variable, there are additional fixed points 
and higher order cycles out in the plane of finite noise amplitude. These fixed and cycle 
points are either hyperbolic or elliptic. In the case of the logistic map we find that each 
standard period-doubling on the x-axis is always associated with a period-doubling out in 
the plane. These cycles in the plane of order 2" are born as elliptic points moving away from 
the cycle points on the x-axis at the bifurcation points. 
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